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Abstract:  We  present  a  B -spline -Bezier  representation  of  r  splines,  curvature  and 
torsion  continuous  quintics  which  have  been  introduced  in  CAGD  bv  Hagen  in 
1985.  Explicit  formulas  are  given  for  the  conversion  of  the  B-spline-Bezier  represen¬ 
tation  to  the  r  spline  representation  and  vice  versa,  and  conditions  and  certain  ranges 
of  tension  values  are  derived  which  insure  the  positivity  of  the  design  parameters. 


0.  Introduction 


In  1974  Nielson  [Nielson  74]  gave  a  piecewise  polynomial  alternative  to  splines  under  tension 
[. Schweikert  66],  [ Cline  74],  the  so-called  v  spline.  The  v  splines  are  curvature  continuous  in¬ 
terpolating  cubics  and  they  are  the  solution  of  the  minimization  of 

ff"||X"(0l|2  dt  +  £v/UX'(',)ll2  V/>0  (1) 


over  the  space 


H 2  =  {  X:  X”  e  L2[t0,  tN],  X' absolutely  continuous  on  [f0,  tN]  } 


subject  to  interpolation  conditions  and  certain  end  conditions.  [Nielson  74,  5<5] 

In  [ Boehm  ?5]  a  B-spline-Bezier  representation  of  v  splines  was  given,  the  y  spline,  and  Boehm 
also  pointed  out  the  close  relation  to  Rarsky's  jl  splines.  [Barsky  81]  Some  details  on  this  relation 
can  be  found  in  [Fritsch  86]  especially  conversion  equations  between  /?  splines  and  v  splines. 
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In  1985  Hagen  [ Hagen  85]  generalized  Nielsons  approach  by  using 


riix<*>(«)u2  dt  +  y,  x  vu  nxa)(f/)ii2  dt  > 

I=0L=\ 


( K  >  2),  rather  than  (1),  minimizing  now  with  respect  to 

HK  =  {X:  X®  e  L2[/o,  t^],  absolutely  continuous  on  [fo>  } 

and  satisfying  interpolation  and  generalized  end  conditions.  Hagens' s  concept  of  'geometric  spline 
curves'  includes  for  K=  2  Nielson's  v  splines  and  yields  for  K  =  3  to  curvature  and  torsion  con¬ 
tinuous  quintics,  the  r  splines. 

The  aim  of  this  paper  is  to  give  a  B-spline-Bezier  representation  of  r  splines.  We  also  derive  the 
conversion  equations  between  the  Bezier  and  the  originally  given  Hermite  representation  in  [// agen 
55].  We  discuss  positivity  conditions  on  the  design  parameters  of  the  Bezier  representation  and 
value  ranges  for  the  point  weights,  the  v,'s. 

Because  we  like  to  fmd  a  B-spline-Bezier  representation  of  r  splines,  we  first  introduce  the  Bezier 
representation  of  segmented  curves.  In  section  II  we  give  a  short  discussion  of  Nielson's  v  splines 
while  Hagens' s  t  splines  are  discussed  in  section  III. 


I.  Bezier  representation  of  segmented  curves 

Let  X(r)  be  a  planar  or  spatial  parametrized  curve  defined  with  respect  to  a  partition  of  the  domain 
space  by  ’knots' 

f0  <  fj  <  ...  tN. 

The  parameter  space  segmentation  induces  a  curve  segmentation  in  Segments 
X,:  [/„  ?K1]  -»■  R*'  (d  =  2,3).  A  local  parameter  u  e  [0,  1]  can  be  introduced  such  that 


/  =  0 . N- 1 


X(0  =  X,(u) 


by  the  linear  interpolation  of  t,  and  tn,: 

t  =  (1  -u)t,  +  utM  , 

The  derivatives  have  to  be  calculated  now  by  the  chain  rule,  i.e. 

=  -trirH*) 

where  A,  =  tM  —  t, . 

Now  the  segments  might  be  given  in  Bezier  representation,  that  means 


for  fe[f„rw] 


where  u  e  [0,  1], 


n 

Xf(w)  =  ^  b„i+kBl(u) 


where  b„/+*  {d  =  2,3)  ,  u  e  [0,1]  and 


=  (£)«*(! -*)"- 


are  the  (ordinary)  Bernstein  polynomials  of  degree  n  in  u.  The  coefficients  by  are  called  Bezier 
points.  They  form  in  their  natural  ordering  given  by  their  subscripts  the  vertices  of  the  so  called 
Bezier  polygon. 
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The  derivatives  of  X,(u)  with  respect  to  u  are  given  by 

irU*  -  g  ^«rci 

where 

A\  =  Ar-I(ba  +  1  —  b  ) 


so  that  m  X(f;)  =  X/_,(l)  =  X,(0)  ,  the  common  boundary  point  of  X/.^u)  and  of  X,(u)  ,  we  have 
for  the  left  sided  derivatives  of  X(f) 


X(r\tJ)  =  lim_  X{r\t)  =  — i 


Arb„ 


A/_,  (n  -  r) !  “ 


and  for  the  right  sided  derivatives 

X(r)(//+)  a  lim  X(rV)  - 


<f)^  =  -1 - A\, . 


r-  n 


A/  («  -  r) ! 


//.  Nu- splines 


Xielsoris  v  splines  are  solutions  of  the  minimization  of  (1)  over  the  space  H2  subject  to  the  in¬ 
terpolation  conditions 


/  =  0 . N 


Xfo)  =  X, 


and  one  of  the  following  end  conditions 

i. )  X'(r0)  =  Xg ,  X'(^)  =  x;, 

ii. )  XM(to)  =  voX'(fo).  X*(r„)  ~  v/vX'(rftr), 

iii. )  X(?0)  =  X(ty) ,  X'(f0)  =  X'l^) ,  X"(f0+)  —  X'(f,v)  —  (vo 


vw)  X'(f0)  • 


=  1, ,  N  -  1)  the  continuity  conditions 

X(tf)  =  X(tf) 

(4) 

x'(tf)  =  xvn 

(5) 

X'(f;+)  =  XVF)  +  v;X'(f/-) 

(6) 

what  can  be  written  in  matrix  form  as 


x;  -  a,x; 


The  (r+  l)2  matrix  At  is  called  connection  matrix  and  the  'vector'  Xf  with 
R*  {d  >  r,  here  is  r  =  2)  is  sometimes  called  the  r-jet  of  X.  . 


n  For 

r+  1  elements  of  4.&I 


■A 


Because  the  curvature  of  a  planar  resp.  spatial  curve  is  given  by 

I  X',  X"| 


ced 

, '  loti _ A 


||X‘ 


»ii3 


resp. 


|[X'  x  X* 
IIX'H3 


(7) 


we  see  that  a  v  spline  is  curvature  continuous. 
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A  Bezier  representation  of  v  splines  can  be  derived  by  inserting  (3)  for  n  =  3  into  (4)  to  (6).  We 
get  as  continuity  conditions  of  the  Bezier  representation  (Figure  1) 

(1  +  <?/)  b7  =  <?/b3/_,  +  b3/+1 


where 


and 


(1  +  y/<7/)b3 /_!  -  '//<?/  b3/_2  +  S; 
(Yi  +  <7/)b3/+1  =  qj  s,  +  y/b3/+2 


4/ 

Yi 


(8) 


/  =  1, ... ,  N  —  1  .  (8)  allows  the  evaluation  of  the  y7 's  of  the  Bezier  representation  of  a  v  spline, 
i.e.  the  evaluation  of  the  y,  's  for  given  v,  values.  On  the  other  side,  the  corresponding  v  spline 
to  a  given  curvature  continuous  cubic  Bezier  spline,  a  so-called  y  spline,  has  v,  values  given  by 


/  =  1 . n  —  1 


(9)  was  first  given  by  Boehm  in  [ Boehm  55].  He  also  presented  a  B-spline  representation  for  cur¬ 
vature  continuous  cubics,  and  pointed  out  the  relation  to  Barsky's  uniformly-shaped  (3  splines. 
[Barsky  5/]  The  connection  to  general  /?  splines,  sometimes  called  explicit,  discrete  or  discretely- 
shaped  (3  splines  [Hoellig  86],  [ Bartels  et  al.  87],  was  also  pointed  out  by  Nielson  and  given  by 
Fritsch  (see  [Fritsch  5(5]). 1  In  fact,  looking  at  the  connection  matrices  of  y  splines,  v  splines  and 
(general)  [3  splines  [Di  n  et  al.  55]  we  see  that  /?,  y  and  v  splines  are  nothing  else  than  different 
representations  of  curvature  continuous  cubics. 


Figure  1.  Construction  of  the  Bezier  polygon  for  v  splines 

Dyn  and  Mkchelli  [ Dyn  et  al.  55]  have  shown  that  the  existence  of  non-negative  local  support  basis 
functions  which  sum  to  one  follows  -  for  geometrically  continuous  spline  curves,  i.e.  spline  curves 
with  continuous  differential  geometric  invariants  like  curvature,  torsion,  etc.  -  from  the  total 
positivity  of  the  connection  matrix.  For  v  splines  the  total  positivity  of  A;  has  the  meaning 
v,  >  0  and  that  is  exactly  the  range  for  the  tension  parameters  v,  covered  by  the  minimum  norm 

1  The  first  local  basis  for  GC1  splines  was  developed  by  Nielson  and  Lewis  in  1975.  [ Lewis  75]  A  local 
basis  for  uniformly -shaped  /?  splines  was  given  in  [ Barsky  81]  and  for  discretely-shaped  /?  splines  in 
[ Bartels  et  al.  84],  see  also  [ Cohen  57], 
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characterization  of  v  splines.  [ Nielson  ;4,  86]  For  v,>0  we  obtain  from  (9)  the  y,  range: 
Yi  <,  1  (remark:  v,  =  0  «=>•  y,  =  1  is  the  usual  C1  cubic  spline).  But  working  with  y  splines  we  know 
that  y,  >  l  are  possible  as  well.  Indeed,  if  we  request  positive  design  parameters,  i.e.  y,  >  0,  so  that 
such  important  properties  like  the  convex  hull  and  the  variation  diminishing  property  are  given,  (S) 
vields  to 


/=  1......V-  1 


V/  >  -  -f"  (1  +  <//)  ■ 

A/ 


Hence  not  only  positive  tension  values  but  v,  values  in  the  range  given  by  (10)  guarantee  positive 
design  parameters  (remark:  v,<  0  yields  to  y,  >  1)  and  therefore  properties  like  the  two  mentioned 
above.  This  result  goes  conform  with  work  done  by  Barsky  [Bar sky  84]  who  extended  the  theory 
of  v  splines  by  identifying  certain  ranges  for  the  v/s  that  guarantee  a  unique  solution  of  the  in¬ 
terpolation  problem.  In  the  special  case  of  a  uniform,  an  equidistant  parametrization,  i.e.  q,  =  1, 
Barsky  gives  the  range  vf  >  —4  which  is  also  given  by  (10). 


III.  Tau-spiines 

Hagens' s  r  splines  are  solutions  of  the  minimization  of  (2)  over  the  space  HK  for  K  =  3  subject 
to  the  interpolation  conditions 

I  =  0, ,  N  X(t,)  -  X7 

and  one  of  the  following  end  conditions  (L  =  K, ... ,  2{K  -  1) ) 

i. )  =  X^-'-"> ,  X<2*->- V)  =  , 

ii. )  X(LVo+)  =  vo^-.-iX^-^Vo) .  x(t)(^)  =  , 

iii. )  X(<b)  =  X(tN) ,  X{2K-l~L){t0)  =  x{2K~l-L\tN) , 

-  X^L\tN)  =  (V0,2K~\~L~  vN,2K-1-D^-  2K  \lo)  ■ 
t  splines  fulfill  at  any  knot  t,  (/=  1, ... ,  /V—  1)  the  continuity  conditions 


x(r/")  =  x(?n 

(ii) 

X'(r;+)  =  XV7) 

(12) 

x"(r,+)  -  x-(fj T) 

(13) 

X'"(f/+)  =  X'VT)  +  *1,2  x'('D 

(14) 

xV)  =  X'VD  -  */.,  X'(0 

(15) 

Because  the  curvature  of  a  spatial  curve  is  given  by  (7)  and  the  torsion  by 

|X',X»,  Xw[ 
l|X'  X  X'll2 

we  see  that  a  r  spline  is  curvature  and  torsion  continuous. 
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The  connection  matrix  of  a  general  torsion  continuous  contact  of  two  curve  segments  of 
S*  (d  >  3)  ol  degree  n  {n  >  4)  -  we  speck  about  geometric  C 3  continuity  .  briefly  GC3  continuity 
-  is  given  by  [ Dyn  el  at.  55],  [ Lasser.Eck  55]: 


Accj 


10  0  0 

0  co  ]  ]  0i  0 

0  a> ;  ^  a>  f  ]  0^ 

0  a>] j 


Comparing  Acci  with  (11)  to  (14).  we  see  that  r  splines  do  not  take  advantage  of  all  shape  pa¬ 
rameters  offered  by  the  concept  of  torsion  continuity.  They  rather  form  a  subset  of  the  set  of  GO 
continuous  curves. 

Furthermore  we  like'to  mention  that  r  splines  and  visual  C3  continuous  curves,  briefly  VC3  con¬ 
tinuous  curves,  i.e.  curs  es  having  contact  of  order  r  [ Geise  62]  with  r  =  3  ,  are  spanning  two  'almost 
totally  separated'  subsets  of  the  set  of  GO  continuous  curves.  The  connection  matrix  of  two  curve 
segments  having  contact  of  order  3  is  given  by  [ Dyn  et  al.  55],  [ Lasser  55]: 

’  1  0  0  o' 

A  0  o,  0  0 

0  u2  o,  03 

0  uj  3u]t>2  i>i 

and  therefore,  comparing  At-Ci  with  A<;Cj,  VO  continuous  curves  form  a  subset  of  the  set  of 
GO  continuous  curves,  and,  comparing  A^cj  with  (1 1)  to  (14),  we  see  that  the  only  curves  beeing 
t  splines  as  well  as  VO  continuous  curves  are  the  usual  O  continuous  curves. 


To  find  a  Bezier  representation  of  r  splines  we  insert  (3)  for  n  =  5  into  (11)  to  (15)  and  get  as 
continuity  conditions  of  the  Bezier  representation  (Figure  2) 


(1  F  qj)  b7  -  q[  b5/_,  F  b5/+1 

(1  F  //<?/)  b5/_j  =  y;< 7,'b5/_2  +  S; 
b  I  +  <7/)b5/+i  =  <7fs/  +  V/b5/+2 

(1  +  Srf,)  b5/_2  =  5tqt  b5/_3  F  e“ 

(<>/  +  £/<7/)s/  =  «/<?/  e7  +  6  fit 
(e7  F  qj)  b5/+2  =  qi*i  F  £/b5/+3 


(i  +  Prq/)^5i-3  =  Piqi^n-i  +  ft 
(pi  F  <7 /<?/)«/  =  ff,q,{J  F  P[t[ 


(<7/  F  r jq,)  ^  = 
(t/  a-  qj)  b5/+3  = 


T,q,t,  F  otf 

+  T/*>5/-Ml 
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*  +  ~  ~ r~  v/.; 

(1-4/)'  3 


4/  A/ 


( 1  +  4/)' 


34/  -  R, 


al  ~  4/'  Z“ 

T,  ~  i  +  (T,-Sj)-±- 


rl  =  4/ 


3  -  r. 


1  4/ 

3<"'  l+  77727 


1 


(1+4/) 


1  -  <7/  A/ 


.  .  24/  A/ 

1  +  - - ttv« 


(1  +  4/) 


3  -  4/  + 


4/  A/ 


(1  +  4/) 


4/  A/ 


( 1  +  4// 


( 20.)  to  (24)  allows  the  evaluation  of  the  design  parameters  of  the  Bezier  representation  of  a  r  spline 
and  by  tliis  the  Bezier  representation  is  given.  The  design  parameters  are  of  course  not  independent 
of  each  other  any  longer.  This  is  obvious  because  two  shape  parameters  are  given  by  (11)  to  (15) 
but  five  are  given  by  (16)  to  (19).  The  dependences  are 


s  _  4/£/ 

C/  1  —  (1  4/)e/ 


resp.  tj  = 


4/  +  (1  “4/)<5/ 
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Figure  2.  Construction  of  the  Bezier  polygon  for  a  r  spline 


(  l  4-  <//0/)( Z;  +  0 t)a i 
3(l4-<?/)<7/  —  ( <3/  +  q  {>-.{) 

{6,  -F  q,zl)T.lTl _ 

1  -r  -  (t/  +  <?/)(£/  +  0;) 

_  ( r-i 

((  l  +  ‘7/<j/)£/ 


_ (<)/  +  q,t[)6,pi _ 

3(  1  +  qi)6jpi  -  (1  +  ‘//<3/)(e/  +  3/) 

1  l£/  +  <?/)(£/  ~h 

£/  3(1  +  <//)<?/  —  (Of  +  <7/£/) 

(I  +  ^0)6/7/ 

(e/  +  <?/)3/ 


That  means  if  an  mterpolating  r  spline  is  given,  i.e.  vM  and  vl2  values,  then  the  design  parameters 
of  the  Bezier  representation  of  the  r  spline  have  to  be  determined  in  such  a  way,  that  they  fulfill  the 
equations  (25)  to  (28)  And  the  equations  (25)  to  (28)  are  valid  if  the  design  parameter  are  deter- 
rnuied  by  (20)  to  (24).  On  the  other  side,  the  design  parameters  of  a  torsion  continuous  quintic 
Bezier  spline  curve  have  to  fulfill  the  conditions  (25)  to  (28)  to  make  the  quintic  Bezier  spline  beeing 
an  mterpolating  r  spline  in  Bezier  representation.  That  means  either  6,  or  z,  can  be  chosen  as 
independent  design  parameter  and.  let's  say  in  case  of  6,,  z ,  has  to  be  determined  by  (25)  as, 
£,.  =  £/(<),)  and  in  addition  either  p,  or  <?,.  or  r,  can  be  chosen  as  independent  design  parameter 
and.  e.g.  in  case  of  p„  i,  and  r,  have  to  be  determined  by  (26)  and  (28)  as  a,  =  a,(p,)  and 
r.  =  T-lpp.  If  the  design  parameters  of  the  Bezier  representation  arc  determined  in  this  way.  then 
the  Bezier  spline  is  equivalent  to  an  interpolating  r  spline  having  the  property  of  minimizing  (2) 
and  we  can  calculate  the  v;  i  and  v;  ;  that  means  the  point  weights  defining  the  jumps  of  the  third 
and  fourth  derivatives  in  the  knots  of  the  r  spline  by 


-7- (I  +<?/)-( -p-  1) 
&r  °i 
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and  so  on. 


The  t  splines  form  a  subset  of  the  set  of  GO  continuous  quintics.  The  B-spline-Bezier  represen¬ 
tation  of  Tau  splines  is  therefore  identical  with  the  B-spline-Bezier  representation  of  GO  contin¬ 
uous  quintics  which  was  given  in  [lick  87],  [Lasser.Eck  §5],  with  the  restriction  that  for  a  t  spline 
the  design  parameters  >■,.  5,.  c,.  pr  a,  and  z,  can  not  be  chosen  independently  of  each  other  as  in 
case  of  the  GO  continuous  qumtic  spline  curve.  We  rather  have  to  set  y,  =  1  and  have  to  choose 
8..  £;.  p-,  o,  and  t,  according  to  the  dependences  given  by  (25)  to  (2S). 


Let's  discuss  now  the  positivity  of  the  design  parameters  for  r  splines.  The  dependences  (25)  to 
(2S)  unply  the  following: 


Because  of  (25),  positivity  of  6,  =  <5 ,(£,)  needs 


0  <  c  <  • 


£  >  0 


in  case  of  q,  <  1  «=•  A,  <  A,_, 
in  case  of  q,  >  1  •=>  A,  >  A,_, 


On  the  other  side,  positivity  of  c,  =  z^o,)  needs 


6  >  0 


0  <  d  < 


in  case  of  q,<  1  «►  A,  <  A,_ 


in  case  of  q,  >  1  »  A,  >  A,_[ 


Because  of  (27)  and  (28),  positivity  of  p,  =  p,(z,)  and  of  o,  =  c,{z,)  needs 

(e/  +  <7/)(£/  +  8j) 

3(1  +  <?/)£/ 

On  the  other  side,  because  of  (26)  and  (27),  positivity  of  p,  —  p,(o,)  and  of  r,  =  t,(<t;)  needs 

8 i  +  <7/C/ 

a  >  — — -  , 

3(  1  +  q/) 

and  because  of  (26)  and  (28),  positivity  of  a ^a^p,)  and  of  r,  =  r ,{p,)  needs 

( 1  +  <7/'5/)(£/  +  Sr) 


H 1  +  <7/1-5/ 
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The  minimum  norm  characterization  of  r  splines  works  for  non-negative  tension  parameter  vn 
and  vt  ,  only.  As  for  v  splines  we  can  extend  the  theory  of  r  splines  by  requesting  the  positivity 
of  the  design  parameters. 

Because  ot  (20)  and  (21).  for  s,  >  0  and  b,  >0  the  tension  parameter  v/2  has  to  be  within  the 
ranee 


that  means 


v/.2  >  max  ( 1  +  4i)2  ~~  ,  -j-1 


v/,2  >  “ 


v/,2  > 


1 


a  +  3 

<7/  A, 
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v/,2  >  ~  ( 1  +  til)1  ~7~ 


s 


if  q,>  l  o  A ,  >  A;_, 
if  <7/  =  1  ~  A,  =  A;_, 
if  q,<\  O  A,  <  A/_, 


Thus  not  only  non-negative  but  also  certain  negative  vu  values  are  allowed. 


(29) 


% 


Because  of  (22)  to  (24),  for  p;>0,  a,>  0  and  r,>0  the  tension  parameter  vf  l  has  to  be  within 
the  ranges  given  by  (30). 


v/.l  <  -fytU  +  <7/)3  +  <7/A,v/2] 
A  / 


v/,t 


>  - 


24  ( 1  +  qtf  +  2(1  -  qj  +  q ’/)  —  vJ  2 


Af 


2  + 


1  +  ^/3 


*/,2 


(30.1) 


(30.2) 


First  v,j  has  to  be  chosen  such  that  ef>0  and  <5;>0  is  fulfilled,  i.e.  v,2  has  to  be  chosen  within 
the  range  given  by  (29),  than  v!t  can  be  chosen  within  the  range  given  above. 

For  q,  =  1  (30)  yields  for  example  to 

<  V;,  <  -^(8  +  A;vf>2) 

A  /  Af 


Barsky  [Bar sky  84]  extended  the  theory  of  v  splines  by  identifying  certain  ranges  for  the  v/s  that 
guarantee  a  unique  solution  of  the  interpolation  problem.  This  idea  allows  especially  the  consider¬ 
ation  of  different  end  conditions.  The  same  can  be  done  for  t  splines,  and  is  indeed  the  topic  of 
actual  research. 

Furthermore  the  idea  of  [Salkauskas  75]  and  [Foley  86,87]  of  introducing  interval  weights  can  be 
picked  up  to  create  interval  weighted  geometric  spline  curves2  minimizing 


^Pi  f  l|X(/O(0H2  dt  + 
/=  0  Vl 


2Zvuiix<%)iia*  • 

/=0  L=  1 


For  A'  =  3  we  get  interval  weighted  r  splines.  Actually  we  are  working  on  this  too. 


2  These  interval  weighted  geometric  spline  curves  are  in  general  not  curvature,  torsion,  etc.  continuous. 
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